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Abstract. We investigate the impact of non-linear corrections on dark energy parameter 
estimation from weak lensing probes. We find that using halofit expressions, suited to 
ACDM models, implies substantial discrepancies with respect to results directly obtained 
from A^-body simulations, when w{z) ^ —1. Discrepancies appear strong when using models 
with w'{z = 0) > 0, as fiducial models; they are however significant even in the neighborhood 
of ACDM, where neglecting the degrees of freedom associated with the DE state equation 
can lead to a misestimate of the matter density parameter i}^. 
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1 Introduction 

The evidence in favor of dark energy (DE), a smooth component with largely negative state 
parameter it; ~ —1, is more than a decade old. The original claim based on the acceleration 
in the expansion rate, deduced from the Hubble diagram of supernovae la [1], has been 
confirmed and strengthened by other probes, in particular cosmic microwave background 
(CMB) spectra (see, e.g., [2]) and large scale structure (LSS) data [3]. The nature of DE is 
however one of the major puzzles of today's physics and astrophysics. In order to investigate 
it, more data constraining the DE state equation are sought. 

Data based on the CMB spectra C^, however, only scarcely constrain the equation of 
state parameter of the Dark Energy component w{z). In fact, although the integrated Sachs 
&: Wolfe (ISW) effect depends on w{z), it mostly affects low~^ which are also affected by 
the ^-dependence of the optical depth r. Furthermore, ISW is an integral effect, and low-^ 
spectral components are subject to a large cosmic variance. 

Putting together CMB and LSS data is already more effective, as this tests the fluctua- 
tion distribution and amplitude at two different z values. Clearly, estimating the fluctuation 
spectrum P{k,z) = {\5Tn{k, z)\^) at several z values, so to work out the growth rate 



would open a real window on w{z) and, therefore, on DE nature. So far G is only very weakly 
constrained (see, e.g., [4]), but a formidable array of surveys are being performed or planned 
to measure it through weak lensing (WL), high-z power spectra and cluster counts [5]. 

This paper will deal with WL probes (see, e.g., [6] for a thorough review). The power 
of this tool is illustrated by the fact that, to fully exploit next generation WL surveys, we 
need predictions on non-linear power spectra accurate up to ~ 1% [7]. 

Such precision goes beyond the claimed ±3 % accuracy of the popular halofit ex- 
pressions [8] based on the halo model of structure formation and calibrated using numerical 
simulations of ACDM models; they yield a map of linear k's onto non-linear /c's, so that 
Pni{k) = Piin{K-)- Leaving apart the need of a generic improvement of halofit accuracy [9], 



G{z,k) = 5^{Kz)/5^{k,Q) 



(1.1) 
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previous work [10] already showed that, using HALOFIT expressions for non-ACDM models, 
requires suitable corrections. In spite of that, the halofit map has been often used to obtain 
the spectra of models with non-constant DE state parameter w{z), that we shall denomi- 
nate dynamical DE (dDE) models, herebelow. This procedure was dictated by the lack of 
appropriate extensions of halofit to non-ACDM cosmologies. 

The aim of this paper is to quantify the effects of using non-linear dDE spectra obtained 
with the HALOFIT map, when the nature of DE is investigated through WL surveys as Euclid 
[11] or COSMOS [12]. 

To this aim we focus on dDE models with 

w{z) = Wo + Wail - a) = Wo + WaZ/{l + z) (1.2) 

(a being the scale factor), although bearing in mind that they are already a restricted class 
of dDE models, their state equations being fully characterized by amplitude and derivative 
at 2; = 0. Through a Fisher Matrix (FM) approach, we evaluate the errors in WL estimates 
for wq, Wa and (the matter density parameter). We do so using matter power spectra 
obtained either from halofit or A^-body simulations and then we compare the likelihood 
ellipses on the parameters obtained with the two approaches. To our knowledge, this is 
the first paper in which a Fisher matrix approach is performed directly on A^-body outputs 
rather than on fits or other approximations. 

We will assume that the shear field is observed in a survey which approximate the 
current design of Euclid [11]. This will be a survey of approximately half sky (20,000 square 
degrees) with galaxies at an average depth of Zm ~ 0.9 and 40 galaxies per square arcminute. 
The redshifts will be evaluated photometrically assuming a normal distribution with variance 
o"2 = 0.05 for errors. The tests are performed assuming as fiducial cosmology ACDM model 
{w = —1), as well as for two dDE models, still consistent with the CMB data from WMAP7 
[2] and related data, roughly the most distant models from ACDM that are still allowed by 
the data (at 95% confidence level (CL) ) in the Wo-Wa plane (see figure 1). In this way we 
explore the dependence of our results on the assumed fiducial model. We can anticipate the 
conclusion that a straightforward use of halofit is significantly misleading, in quite a few 
cases. 

2 Weak lensing parameter estimation 

The first step to estimate the errors in the parameter measured through weak lensing tech- 
niques, amounts to define the convergence weak lensing power spectrum, which in the linear 
regime is identical to the ellipticity power spectrum [6]. This spectrum is a linear function 
of the matter power spectrum convolved with the lensing properties of space. According to 
[13], if we bin the galaxies in rif, redshift bins labeled by i,j = 1, ...n;,, it reads 

Here Pni{k, z) is the non-linear matter power spectrum at redshift z and the Wi are window 
functions that will be explicitely given below. Through this paper rib = 3 and 5 will be 
considered, in order to assess the dependence of the results on the binning procedure (there is 
no much gain in going beyond 5 bins [14]). In the last two cases, the bin limits Zi are selected 
so to have the same number of galaxies per bin. Let us then define the other quantities in this 
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Figure 1. Constraints on tlie time-dependent DE equation of state (1.2), derived from WMAP7 
spectra combined with SNIa and BAO updated data [2]. The contours show the 68% and 95% 
CL. The square point indicate the ACDM model; the triangle points are the other two reference 
cosmologies considered in this work: the model with negative (positive) Wa is also dubbed Ml (M3). 
The intersection of the dotted lines is the best fit point in [2] Wo = —0.93, Wa = —0.41. 

expression, starting from the redshift dependence of the Hubble parameter H{z) = E{z)Hq 
and 

accordingly, by referring to the scale factor a and the conformal time r, we have that 

dzlE{z) = HQ{-da/a^)/{h/a^) = -Hodr , (2.3) 
so that eq. (2.2) also reads 

r{zi,Z2) = Ho{n-T2) . (2.4) 

From here we deduce that r{z) = r(0, z) = HqItq — t{z)]. 

Let us then consider also the distribution of the galaxy number per unit redshift and 
solid angle. We assume the following simple parametrized form 

n{z) = d^N/dn dz = C {z/zo)^exp[-{z/zo f] with C = (B/zq) r-^[{A + l)/B] (2.5) 

where T{x) is the Gamma function. Here we choose A = 2, B = 1.5, so that C = 1.5/zo. If 
Zm = 0.9 is the average redshift, one then has zq = Zm/1-412. This choice of values is a good 
approximation to the observed selection function of a survey like Euclid [11]. By integrating 
the distribution n{z) across the depth of the j-th bin Azj we have the angular density in 
steradians 

Uj = 3600(180/7r)2 / dz n{z) , (2.6) 

JAzj 
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Figure 2. Di{z) distributions in the 3-bin case. The bin limits zi & 22, shown in the frame, are 
selected to have the same number of galaxies in each bin. In the 5-bin case the bin limits zi are 0.560, 
0.789, 1.019, 1.324, 3.00 . The background black line is n{z). 

an expression to be used in the next section. 

The distribution (2.5) is then considered within the limits of the redshift bins. In order 
to take into account the discrepancies between the photometric redshift in use, and the actual 
galaxy redshift, we filter the distributions using window functions [15, 16] 



with (t{z) = Uz (1 + -z), o"^ = 0.05 (see, e.g., [17, 18], for the motivation of this parameter 
choice) . 

In figure 2 we show the radial galaxy distribution and the distributions 



in the case of three bins. They must then be normalized, yielding the distributions 





(2.7) 
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wherefrom we derive the functions 




(2.10) 



shown in figure 3. 
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Figure 3. Wi{z) window functions adopted for the 3-bin case. 

From the power spectrum Pij{£), we obtain the covariance matrix 

Cjk = Pjk + Sjk{llt)nj' , (2.11) 

1 /2 

including the effect of the r.m.s. intrinsic shear jint, for which we assume {^fnt) = 0-22 
[17]; Uj is given by eq. (2.6). 

The above expressions enable us to make use of the Fisher matrix formalism, which 
provides lower limits to the errors on the cosmological parameters one aims to measure. Its 
basic tool is the likelihood function, yielding the probability that a model, fixed by a set of 
parameters pa, gives the set of data x. 

The Fisher matrix for tomographic weak lensing reads [13] 

- fsky ^ ^ ^C^-fc -^Cmi ■ (2.12) 

The basic ingredients of this expressions are the Cij components; fs^y is the sky fraction 
covered by the experiment (in our case f^^^y = 0.5); the cosmological parameters pa we shall 
consider are Wo, Wa and Q,m, so that d/dpa are partial derivatives with respect to them. 

All the other parameters are fixed, in order to keep the number of A^-body simulations 
to a manageable level, by assuming data to provide us the other parameters with negligible 
errors. These parameters include O;,, n^, h (baryon density parameter, primeval spectral 
index of scalar fluctuations, adimensional Hubble parameter), plus a parameter setting the 
spectral normalization. 

If data yield non-linearly evolved amplitudes, it is from them that we should guess 
suitable initial conditions at = 24, when starting the simulations, in the linear regime. 
Here we shall need to improve on the standard normalization procedure which neglects the 
mild discrepancy between the r.m.s. matter density fluctuation on the comoving scale R% = 
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8 h Mpc, taken from data, and the linear expression 

POO 

= (27r2)-W dk Ak''^+^T\k,z)W\kR8) ; (2.13) 
Jo 

here A is the primeval fluctuation amplitude at a suitable very high zq value, T{k, z) is the 
linear transfer function from such zq value to z, W{x) is a window function, for which we 
can use the Fourier transform of a top-hat filter: W{x) = (3/x^)(sinx — xcosx) . This mild 
discrepancy, infact, is 0(10%), too much for the signal we are chasing. 

Normalizing models to have the same cr8,n/) as we then need to do, is clearly harder and 
includes technical aspects which will be further discussed below. 

Should we forget this point, and normalize all models so to have the same erg,; at z = 0, 
the derivatives with respect to Wq, Wa or Vtm are largely dominated by the normalization shift 
at ^ = 0, as the a^^i-a^^ni shift itself depends on Wq-, Wa and Vim- This would confuse the z 
dependence of the growth factor, through the observational z-range, i.e. the main observable 
that lensing measures want to exploit. 

Once the fiducial models are chosen, the derivatives in eq. (2.12) are evaluated by ex- 
tracting the power spectra from the simulations of models close to the fiducial ones, obtained 
by considering parameter increments ±5 %. 

Our task will then amount to find the regions above a given confidence level. They are 
M-dimensional ellipsoids, about the M parameter values corresponding to fiducial models, 
that we shall then project onto 2-dimensional subspaces spanned by parameter pairs, so 
determining constant contours. In the Fisher Matrix that we adopt, such marginalized 
likelihood contours will be ellipses centered on the fiducial cosmology. 

In a generic case, when the M parameters are dubbed {a = 1,2, ..,M), the ellipse 
equation in the , plane (n = 2 subspace) reads 

(^£^=^x\CL,n) . (2.14) 

Here Ax^ sets the confidence level, being 

Ax^(68.3%, Up) = 2.30 , Ax^(95.4%, rip) = 4.61 , (2.15) 

for Up = 2 . 

3 Cosmological models 

In figure 1 the settings of the fiducial models are shown. The corresponding values of Wq 
and Wa, are reported in Table 1. We also take 0^ = 0.274, h = 0.7, as,i = 0.81, n = 0.96 
(matter density parameter, Hubble parameter in units of 100 km/s/Mpc, linear r.m.s matter 
fiuctuation amplitude at z = on the scale of 8 h~^Mpc, primeval spectral index, respec- 
tively) consistently with WMAP7 outputs. We shall discuss below why (slightly) different 
as^i values need also be considered. 

About the fiducial models we select models characterized by ±5 % increments of Wq and 
Wa, to perform numerical derivatives. In Table 1 the list of these models is reported. 

When ACDM is the fiducial model, we consider also derivatives with respect to the 
parameter. In this case all the procedure deals with ACDM cosmologies; we shall then work 
out the spectra also for ACDM models with Qm = 0.2603 and 0.2877. 



(XpX,) 
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Table 1. List of values of Wo and Wa used to perform numerical derivatives. 
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Figure 4. i?-dependence of the ratio between the (T^(i?)'s at z = obtained from simulations, for 
models in the neighborhood of Ml; they arc normalized to the respective mass variances obtained 
using HALOFIT, however almost identical for all these models (see text). 

For ACDM and all nearby models, two different initial seeds were also considered, to test 
the dependence on initial conditions, finding that Fisher Matrix results are almost insensitive 
to it. Accordingly, for the other fiducial models, only one seed is used. Altogether we run 24 
model simulations. 

4 Simulations and their analysis 

A^-body simulations are performed by using a modified version of PKDGRAV [19] able to 
handle any DE state equation w{a). To this aim, A^^ = 256^ particles, representing CDM 
and baryons, are set in a box with side L = 256/i~^Mpc. Particle masses are therefore nic = 
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Figure 5. Connection between halofit and simulation spectra at the low-fc end. The upper (lower) 
panel shows the resulting spectra when normalizing to <7s,i{z — 0) {crs^niiz = 0)). The dotted line is 

HALOFIT. 



Po,cr^m{L/N)^ = 7.61 X 10^^h~^MQ {po,cr '■ Critical density). Softening is 1/40 of the initial 
intra-particle distance, yielding e ~ 25 h~^kpc (wavenumber k = 2-K/e ~ 150/iMpc~^). 

Transfer functions generated using the CAMB package are employed to create initial 
conditions, with a modified version of the PM software by Klypin and Holzmann [20], also 
able to handle suitable parameterizations of DE. Simulations are started at Zin = 24. The 
initial density field is obtained from the same random numbers for all the simulations. For 
the ACDM fiducial we generate each simulation with two sets of initial random numbers. 

Matter power spectra are obtained by performing a FFT (Fast Fourier Transform) of the 
matter density fields, that we computed from the particles distribution, through a Cloud-in- 
Cell algorithm, by using a regular grid with Ng = 2048. This allows us to obtain non-linear 
spectra in a large /c -interval. 

In particular, our resolution allows to work out spectra up to A; ~ 10 /iMpc~^. However, 
for k > 2-3/iMpc~^ neglecting baryon physics is no longer accurate [21]. Although both 
simulations and halofit refer to purely gravitational dynamics, we wish to examine only 
physically significant scales. Therefore, we consider WL spectra for I < 2000 only, so limiting 
the contribution of scales k > 2/iMpc~^ below ~ 5% , at most. 

Let us now discuss the normalization problem. In Figure 4 we show the dependence on 
scale of the r.m.s. of matter density fluctuations (mass variance), as obtainable through the 
expression 

a'^{R) = "^AkkP{k)k^W^{kR) (4.1) 

k 

where P{k) is the full non-linear spectrum. When applied to simulations, we denote the 
variance as a1^^{R). Plotting the scale dependence for a1^^{R)^ rather than for the simulation 
spectra Pgimik), smears out the numerical oscillations that spectra unavoidably exhibit. 
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Figure 4 allows us to focus on the discrepances between the model Ml and the nearby 
models used to perform derivatives about it. The spectra shown in the figure were obtained by 
normalizing initial conditions so that their ag^i at z = coincide. As a consequence the whole 
linear spectra of Ml and nearby models (not shown in the Figure) are nearly coincident, at 
z = 0, since the models share all parameters except the DE state equation, which has a very 
weak impact on the linear spectrum shape. If one uses halofit to correct the linear spectra, 
the z = non-linear spectra would still coincide, since halofit only depends on a^^i and on 
Orfe, not on w. In contrast. Figure 4 confirms that the simulation spectra differ. Accordingly, 
taking the Fisher matrix derivatives fixing ug^; is different from fixing ag^ni- Besides of the 
points made above, here we choose to normalize to the same non-linear variance erg also 
because this is the primary observational quantity and we expect that future data will put 
direct constraints on it, resulting in a strong prior in the Fisher matrix. 

The price to pay is that producing simulations with different cosmological parameters 
but coincident crg^nl is not a trivial task and requires a trial-and-error search of initial condi- 
tions for the A^-body simulations. 

There is then another question, concerning the low-A; spectrum. The point is that the 
integration (2.1) unavoidably extends down to /c's which are still close to (or fully in) the 
linear regime. In principle, dealing with such scales should not create difficulties. However, 
the spectra obtained from our simulations, run in a box of 256 /i^^Mpc, become increasingly 
noisy and discrete when the linear regime is approached. To solve this problem we extrapolate 
our numerical spectra with the linear theoretical ones at scales larger than a certain threshold. 

The problem is then where and how to connect numerical spectra with theoretical ones. 
Let then be kp = 2p7r/(256 /i~^Mpc), being p an integer; we use simulation spectra when 
Akp/kp shifts below 0.1 , i.e. for log(A;//iMpc~^) > —0.61 . At lower /c's halofit spectra 
are used. 

Figure 5 (upper frame) is the result of this operation, when done using M3 ± Awa 
spectra characterized by the same crs,iiz = 0). Spectral shifts of this kind affect the ratios 
between derivatives with respect to cosmological parameters (in this case, Wa)- In the same 
Figure (lower frame) we show the result of the connection after renormalizing spectra to the 
same as^ni- Some scale dependent discrepancy still remains, but the main signal comes then 
from the different redshift dependence of model spectra, rather than from the shift among 

5 Predictions on parameter errors 
5.1 Derivatives 

The next step is evaluating the derivatives in eq. (2.12). For the sake of illustration we plot 
the results of such differentiation in the case of a single bin, by considering: (i) the spectra 
obtained from simulations, as illustrated in previous section; (ii) HALOFIT spectra; (iii) purely 
linear spectra. 

In figure 6 we show the derivatives and, in the lower frame, compare the results of the 
cases (i) and (ii), when the fiducial model is ACDM. The results obtained starting simulations 
from two different seeds are also compared. In figure 7 we similarly show derivatives when 
the fiducial model is Ml or M3. One seed only is considered there, and we consider the 
parameters wq and Wa only. 

The effects of BAO's as well of irregularities in simulation spectra are clearly visible 
in all plots. There is however a clear difference between the derivative with respect to ^Im, 
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Figure 6. Derivatives of C^, in the case of a single bin, for 20 < £ < 2000, with respect to the 
parameters flm, wq and Wa, when the fiducial model is ACDM. In the upper frame we show their 
different behavior, for A^-body simulations, halofit and linear spectra. In the lower frame we show 
the ratio between derivatives obtained from A^-body simulations and halofit. For these models we 
run simulations with two different seeds. Derivatives obtained from each of them are in the same 
colors, using a darker tonality for the seeds considered just for these models. 

in figure 6, and any other one. For both seeds the ratio between its value for simulations 
and HALOFIT is close to unity, with a discrepancy exceeding ~ 5 % just for i > 1500, when 
contributions from k values above the range of HALOFIT validity become non-irrelevant. 
This is a further confirmation of the coherence between A^~body simulations and halofit, 
for ACDM models, the only ones involved in the differentiation. 

Still in figure 6, we notice that differences between seeds mostly stays within 6-7%, 
even in the case of the parameter Wa, which is the most sensitive to noise effects. 

The case Ml, shown in figure 7, exhibits then a peculiarity, with respect to all other 
cases. Here, the ratios of derivatives with respect to wq and Wa almost overlap up to ^ ~ 500, 
while their trend remains similar even up to £ = 1000 . The effect is peculiar of Ml, and 
recurs also when more bins are considered. For ACDM, the behaviors of the two ratios 
are significantly discrepant. In the case M3, finally, the two ratios exhibit neatly different 
behaviors. This reflects onto the shape of the likelihood ellipses which are one of the main 
results of our work. 
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Figure 7. Derivatives of Ci, as in the previous figure, for Ml and M3 fiducial models (l.h.s. & r.h.s., 
respectively), omitting the parameter $7^, which is not considered in these cases. 



5.2 Likelihood ellipses: ACDM neighborhood 

In figures 8 and 9 we show the confidence ellipses, when the fiducial model is ACDM, in the 
cases of 3 or 5 bins and using up to ^ = 500, 1000 or 2000. 

Up to £ ~ 500, i.e. including just a mildly non-linear scale range, discrepancies between 
HALOFIT and simulations already indicate an underestimate of errors in the halofit case. 
Discrepancies approach a factor 2 when we go up to £ = 1000, both for 3 and 5 bins. They 
become quite significant in the i = 2000 case. In the latter cases we observe a modification 
of the likelihood contours, indicating a different correlation between the errors on different 
parameters. This includes the Wa-^m correlation and means that a neglect of the degrees of 
freedom of DE state equation leads to errors on the estimate of the density parameter r^^. 

Meanwhile, the discrepancy between simulations started from different seeds are much 
smaller. Different seeds yield discrepancies still visible in derivatives (see figure 6), but 
the integrations needed to pass from them to likelihood ellipses smear them out. In fact, 
the shaping of the ellipses is dominated by the 2;-dependence of the spectral behavior, and 
here is where halofit is apparently unable to meet simulation results. The changes in the 
confidence ellipses are also quantified in Table 2. 

Incidentally, let us notice that the improvements, when passing from the 3- to the 5-bin 
case, are more significant when a greater I is considered. 

In particular, in Table 2, we give the inverse area of the ellipses at 95% CL ; when 
considered in the WQ-Wa plane, it coincides with the figure of merit 
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3.65 xlO"^ 
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6.57 xlO-2 


9.36 xlO-2 






1.21 xlQ-'-^ 
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Table 2. Features of likelihood ellipses in the neighborhood of ACDM, for the 5 bin case with 
£ = 2000; the angle between the major axis and the abscissa is dubbed 9; the abscissa being wq 
(former two cases) or Wa- For the definition of the inverse area and correlation see text. 

introduced in the report of the Dark Energy Task Force [22]; the factor 4.61 yields 95% CL 
(see eq. (2.15)). 

The correlation coefficient or Pearson correlation measures the correlation between two 
variables X, Y. In terms of Fisher Matrix components, it reads 

p{X, Y) = F^'y/^Jp^l^xFy^ . (5.2) 

Table 2 confirms that, when considering models different from ACDM, non linear cor- 
rection obtained through halofit may be misleading. This is true even when the fiducial 
model is ACDM itself and we just consider mild deviations of w from —1. Let us point out, 
in particular, the Wa-^^m case; here the parameter correlation increases by a factor ~ 5, when 
passing from HALOFIT to simulations, because of a radical change in the orientation of the 
ellipse, as the angle between its major axis and the Wa axis passes from -0.15° to -0.49° (see 
also figure 9). 

As expected, the error on estimate is not affected by the passage from simulations 
to HALOFIT; in this case we are dealing with ACDM models only. On the contrary, using 
HALOFIT leads to underestimates of the errors on wq and Wa, by a substantial 30-40% . 
Let us however outline that the substantial shift of the O^-w^a correlation means that the 
very estimate of the matter density parameter 

is in jeopardy, as soon as we admit mild 
deviations from pure ACDM, so that the density of DE has a slight dependence on redshift. 

5.3 Confidence ellipses: Ml & M3 models 

Figures 10 and 11 then show the results in the WQ-Wa plane, when the fiducial models are 
Ml or M3. It is evident that the two cases are quite different. This refiects the behaviors of 
derivatives, whose ratios almost overlap in the Ml case, while remaining different in the M3 
case. 
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Figure 8. Likelihood ellipses, showing 65 % and 95 % CL contours, for all parameter pairs, including 
signals up to either £ ~ 500 or 1000 for the ACDM fiducial. In the former case no significant 
discrepancy between halofit and simulation results can be appreciated. Discrepancies become more 
relevant in the latter case, both for 3 and 5 bins. Notice also that the two simulation seeds yield 
almost overlapping outputs. 

In the Ml case, we see just quite a mild shift. For i = 500 or 1000, and 3 or 5 bins, 
the discrepancy appears rather small, being just slightly more evident for i = 2000. Let us 
however outline that such mild discrepancies are O (10%), compatible with discrepancies 
attaining ~ 20 % in the derivatives (see figure 7, left panel) . It is also clear that including 
larger ts means integrating over slightly greater discrepancies. The key issue, however, is 
that, until the two derivative ratios proceed with similar trends, the only effect on ellipses is 
a change in their area, while the orientation of their axes does not change. 

Let us now come to the M3 case, where the derivative ratios exhibit radically differ- 
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Figure 9. Likelihood ellipses, showing 65% and 95% CL contours, for all parameter pairs, consid- 
ered including signals up to £ ~ 2000 for the ACDM fiducial. Here simulations and halofit yield 
significantly different outputs. 

ent trends. The most immediate effect is that errors estimated tlirougli halofit exceed 
simulation errors by a substantial factor. Discrepancies are already visible when only mildly 
non-linear contributions are included and increase sharply when exploiting more deeply non- 
linear spectral areas. 
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In Table 3 we report quantitative information on the deformation of likelihood elhpses 
on the 5-bin i = 2000 case. The only elements which do not change much are the correlation 
and, accordingly, the angle 9. On the contrary, the ratio between the axes of the ellipses and 



- 15 - 







HALOFIT 


sim. average 


Wo-Wa 


e 


-80.6" 


-83.9° 


plane 


axial ratio (95% CL) 


49.7 


18.3 




l/area(95% CL) 


181 


505 




correlation 


-0.99 
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0.136 
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Table 3. Features of likelihood ellipses when the fiducial model is M3, for the 5 bin case with i = 2000; 
symbols as in Table 2. 

the FOM vary by a factor ~ 2.5 . Standard deviations are also deeply affected. 

Altogether, this is a case when estimates based on halofit are simply not trustworthy. 

6 Conclusions 

In the last decades, cosmological observations led to the most important discoveries in fun- 
damental physics. The dark cosmic components are either a proof of physics beyond the 
standard model of elementary interactions or a way to parametrize violations of General 
Relativity. In the former case, the two dark components could be the phenomenological 
description of a single cosmic fluid [23]. Energy exchanges between such components could 
then be a signal of such dark unity [24]. In turn, this opens a list of possibilities which is 
embarrassingly long. 

This florilegium of options is a clear consequence of the lack of data. Tomographic 
weak lensing measures will clearly contribute to provide more of them. This undoubtedly 
implies an important experimental effort and there are scarce doubts that, accordingly, data 
interpretation tools deserve to be updated. 

In this work we outline the need to improve our capacity to predict the spectra of 
cosmological models. When we restrict ourselves within the frame of ACDM cosmologies, 
e.g. to perform derivatives in respect to the halofit expressions still yield acceptable 
outputs (see Figure 6). Improving their accuracy keeps however a need, namely for the 
parameter range consistent with most recent observations. 

A first technical conclusion is that an inspection of the non-linear part of spectra, 
through tomographic weak lensing surveys, greatly improves their discriminatory power. 
This is also confirmed by the effects of including an increasing number of i harmonics in the 
analysis: when halofit and simulations yield discrepant results, the difference between them 
increases when including greater and greater i values. On the contrary, with the instrumental 
features assumed here, essentially corresponding to those of the Euclid [11] project, there 
seems to be no clear advantage in going from 3 to 5 bins. 

Our main conclusion is however that, going beyond a DE state equation with w = — 1, 
direct model simulation remains the only efficient way to obtain spectra. Some attempt to 
generalize halofit to constant w ^ —1, although important, do not hit the heart of the 
problem, as we rather expect a variable w{z). We should also be aware of the severe bias on 
w{z) detection which could follow the assumption of constant w [25]. 

In this work we have tested the effects of using ACDM-suited halofit expressions, 
when the true cosmology is a simple dDE model, straightforwardly consistent with WMAP7 
data [2]. 
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Figure 11. Contours at 68% and 95% CL for the M3 model. 



We find that, even with ACDM as fiducial, when trying to estimate the error on wq 
and Wa, we could be misled, if using halofit expressions. Not only they lead there to error 
misestimates, but they can also imply a wrong correlation between the estimate of and, 
e.g., Wa- Henceforth, the very Qm estimate is in jeopardy, if we do not treat accurately 
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the functional ensemble of state equations w{z) about w = —1. This is, perhaps, the most 
striking result of this work. 

We extended our test to the two dynamical DE models, differing from ACDM just 
for w{z), which are however within 95% CL from ACDM, on the basis of available data. 
Curiously enough, the two models behave rather differently. 

For one of them, the discrepancies between halofit and simulations, although visible 
in the derivatives, leave just O (10%) discrepancies in similarly oriented likelihood ellipses. 
In this case error estimates performed through halofit turn out to be fair. 

In the case of the other fiducial model, we meet an absolutely different situation. First 
of all, using halofit leads to error overestimates by large factors (100 to 200% larger), to 
a wrong figure of merit, and to other problems which can be summarized by stating that 
HALOFiT-based results are wrong. 

In the latter cases (non-ACDM) we have not yet tested the intercorrelation between 
DE state equations and the density parameter. This will be among our future tasks. 

In spite of that, these results suggest a trend in the effectiveness of HALOFIT, when non- 
ACDM cosmologies are involved. We meet: (i) An apparent reliability of such expressions, 
when the fiducial model is characterized by negative Wa- (ii) Significant difficulties for the 
ACDM case {wa = 0), where we expected them to have their best performance, (iii) In the 
case of positive Wa, finally, we conclude that halofit expressions yield misleading results. 

We might then conjecture that an indicator of halofit efficiency is the sign of the 
derivative 

w'{0) = ^ (z = 0) . (6.1) 

When w'{0) > 0, any use of HALOFIT expressions seems dangerous. However, as we have 
been dealing only with DE state equations which can be fully characterized at z = 0, this 
conjecture need to be tested for other cases. 

Altogether, let us conclude that halofit, in its present form, is not adequate to deal 
with the degrees of freedom opening when non-ACDM cosmologies are considered, as none 
can predict, a priori, which can be the physical w{z) behavior. 
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